Auxiliary-Fermion Approach to Critical Fluctuations 
in the 2D Quantum AF Heisenberg Model 
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The nearest-neighbor quantum-antiferromagnetic (AF) Heisenberg model for spin 1/2 on a two- 
dimensional square lattice is studied in the auxiliary-fermion representation. Expressing spin oper- 
ators by canonical fermionic particles requires a constraint on the fermion charge Qi = 1 on each 
lattice site i , which is imposed approximately through the thermal average. The resulting interacting 
fermion system is first treated in mean-field theory (MFT) , which yields an AF ordered ground state 
and spin waves in quantitative agreement with conventional spin- wave theory. At finite temperature 
a self-consistent approximation beyond mean field is required in order to fulfill the Mermin- Wagner 
theorem. We first discuss a fully self-consistent approximation, where fermions are renormalized due 
to fluctuations of their spin density, in close analogy to FLEX. While static properties like the corre- 
lation length, (,{T) oc exp(a J/T) , come out correctly, the dynamical response lacks the magnon-like 
peaks which would reflect the appearance of short-range order at low T . This drawback, which is 
caused by overdamping, is overcome in a 'minimal self-consistent approximation' (MSG A), which 
we derive from the equations of motion. The MSCA features dynamical scaling at small energy and 
temperature and is qualitatively correct both in the regime of order-parameter relaxation at long 
wavelengths A > ^ and in the short-range-order regime at A < ^ . We also discuss the impact of 
vertex corrections and the problem of pseudo-gap formation in the single-particle density of states 
due to long-range fluctuations. Finally we show that the (short-range) magnetic order in MFT and 
MSCA helps to fulfill the constraint on the local fermion occupancy. 
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I. INTRODUCTION 



Antiferromagnetic correlations play an important role 
in the physics of high-temperature superconductors. 
Magnetism is most pronounced in the undoped par- 
ent compounds, where the almost uncoupled two- 
dimensional CuO-planes form magnetic Mott insulators 
with one spin-1/2 electron per Cu-site. The simplest 
model for a CuO-plane is the quantum Heisenberg anti- 
ferromagnet (QHAF) in two dimensions (2D) on a square 
lattice. 
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The super-exchange coupling = J > is non-zero 
only if the Cu-sites i,j are nearest neighbors. In the 
undoped compounds studied experimentally, it assumes 
values around J ~ 130 meV. 

This model has frequently been considered in the pas1*i. 
Its ground state at zero temperature T = appears to 
be a Neel state with an order parameter |(S)| moder- 
ately reduced from its classical value 1/2 by quantum 
fluctuations. At finite temperature T > the theorem 
of Mermin and Wagner prohibits any long-range order. 
For sufficiently low T the magnetic correlation length ^ 
grows exponentially in the two-dimensional system and 
diverges at T = , 



^(T) cx exp(a J/T) 



(2) 



^ is measured in units of the lattice spacing, which is 
set to 1 . In the following we also use h = 1 , fc^ = 



1 . The exponential factor in £,{T) has been obtained 
from quantum-monte-carlo (QMC) calculationSi^ of the 
QHAF, high-temperature series expansion^, Schwinger- 
boson mean- field theorji& (SBMFT), and modified spin- 
wave theorj^. The exponent comes ouljSi^ as a ~ 1.15 . 
An exponential behavior like Eq. ^ also results from the 
(classical) non-linear sigma modelM (NLcrM) as well as 
the quantum NLcrM (QNLcrM) in the so-called renormal- 
ized classical regimeiSiiiii^. 

The theoretical predictions from the QHAF and the 
(Q)NL(tM are consistent with experiments on undoped 
cuprates above their respective Neel temperature Tn , 
confirming the almost two-dimensional nature of fluc- 
tuations in these systems. In particular, a correlation 
length behaving like Eq. Q has been seen in quasi-elastic 
(ener gy i ntegrated) neutron scattering on Sr2Cu304Cl2 
iRef. M), SraCuOsCls (Ref. El, and La2Cu04 (Ref. 
Il5lllari7|) . The observedi^ti^ii^ exponent a = 1.15 agrees 
with theory. A three-dimensional (algebraic) behavior 
might be observablei^ii^ at temperatures very close to 
Tn- 

The dynamics of spin fluctuations are measured 
through the dynamical structure factor S{q, to) . In gen- 
eral two regimes are expected, depending on the wave 
vector k = q — Q relative to the Neel ordering vector 
Q = (7r,7r) : For /c^ <C 1 , i.e., for distances much longer 
than the correlation length, the spectrum S shows a sin- 
gle peak at cij = with some width ^ luq that depends 
weakly on A;^ . The energy scale loq represents the time 
scale for relaxation of the order parameter in the dis- 
ordered phase T > . Critical slowing down leads to 
a vanishing of luq like wo~C~^asT^O. The dy- 
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namical exponent z comes out as z = 1 . For fc^ ^ 1 the 
peak position shifts to finite w and becomes /c-dependent, 
reminiscent of damped spin waves. These originate from 
the short-range magnetic order, visible only at distances 
small compared to the correlation length. Moreover, the 
dynamical scaling hypothesisii states that for ^ S> 1 and 
all /c <C TT spin fluctuations are governed by the single 
energy scale ; the structure factor obeys the scaling 
form 

5(q,L.) = — 5^*(q)<i>(fce,c./wo) (3) 

with the static structure factor (equal-time correlation 
function) 

5^*(q) = (5^5^q) =5^*(QMfcO- (4) 

Lp and $ denote (a priori unknown) scaling functions. 
The dynamical structure factor has been calculated^^ 
using the QMC plus maximum-entropy method for 
the QHAF, molecular dynamics (MD) on the classical 
Heisenberg model^^, the SBMFT^O-^i, y^V) on a clas- 
sical lattice modelSS equivalent to the QNLctM in the 
renormalized classical regimeii. The form (PJ of the 
structure factor, the presence of propagating modes in 
the 'shape function' <I> for fe^ > 1 , and the energy scale 
'^o ~ have been well confirmed. The latter has 
also been calculated from coupled-mode theoryS^. Ex- 
perimentally, inelastic neutron scattering (INS) on sev- 
eral cuprate parent compounds--'^^^ at temperatures 
T > Tat shows a quasi-elastic peak at zero energy, the 
observed-^ dynamical exponent is z = 1 . The latter is 
also found in magnetic resonance^i. The INS data is con- 
sistent to spin-wave like excitations at higher energies24 
in the paramagnetic phase. 

Results for the dynamics of spin fluctuations come 
primarily from purely numerical approaches like QMC 
and MD . For an analytical treatment of the Heisenberg 
model, the non-canonical commutation relations of spin 
operators 8^,3^,3^ pose a severe difficulty. This can 
be circumvented by re-writing S'^ in canonical boson or 
fermion creation and annihilation operators. These act 
on an enlarged Hilbert space that contains unphysical 
states, which have to be removed by imposing a con- 
straint. Using Schwinger bosons 6„ , n — 1,2, the spin- 
1/2 operator and constraint read& 

n,n' = l,2 n=l,2 

Alternatively, one can use Abrikosov auxiliary 
fermions^*'^^ /ti/i. that carry spin-1/2, 

S.-^ E , Q.=E = ! (5) 

a,Qi'— ±1 a 

Here a lattice-site index i has been added; a = +1, —1 
corresponds to "f , | . 



The fermion representation |(SJ of the Heisenberg 
model CP is actually the 1/2-filled (undoped) limit of the 
standard slave-boson formulation of the J-model. The 
latter has been extensively studied using, e.g., mean-field 
theory''*'. These works focus on the superconducting and 
normal-state properties of the t-J-model at finite hole 
filling (doping), and its application to the doped cuprate 
superconductor^^. At low hole concentration and in par- 
ticular at 1/2-filling the slave-boson mean-field theory 
gives antiferromagnetic (AF) long-range orde]4i*2S. The 
strong magnetic correlations at finite temperature in the 
2D system, however, have not been considered in this 
type of approach. 

The purpose of this paper is to study the antiferromag- 
netic spin-1/2 Heisenberg model (QHAF) on a strictly 
two dimensional square lattice at finite temperature, us- 
ing the auxiliary-fermion representation (j^l . In con- 
trast to the boson approachSs^, mean-field theory with 
fermions leads to a finite Neel temperature Tat = 0.5 J. 
While mean-field results are good in the AF ordered 
ground state (see Sect. IIII Al below) . at T > a self- 
consistent T-matrix approximation beyond mean-field is 
required in order to suppress the phase transition. Two 
approximation schemes will be proposed and compared 
to the established results reviewed above. A suitable ap- 
proximation should be able to fulfill the Mermin- Wagner 
theorem and reproduce the exponentially growing mag- 
netic correlation length, also it should yield a qualita- 
tively correct description of the dynamics including crit- 
ical slowing down and short-range order. Once the im- 
portant elements of such an approximation are identified, 
the theory can possibly be extended to less well under- 
stood cases, e.g., the i-J-model at finite hole filling by 
introducing 'slave' bosons for the doped holes. 

The present study is related to certain self-consistent 
diagrammatic schemes for the 2D Hubbard model, like 
the FLEX approximation^^: When the magnetic corre- 
lation length is growing (e.g., due to a reduction of hole 
filling) the density of states of electrons is expected to de- 
velop a pseudo precursor gap^* around the Fermi level, 
reflecting the proximity to an AF ordered state which is 
characterized by zero spectral weight at the Fermi en- 
ergy. It has been demonstrated^ that an approxima- 
tion with 'reduced self-consistency', i.e., where part of 
the electron lines in the self energy are replaced by bare 
ones, favor the formation of such a spectral-weight sup- 
pression. A similar observation is made here: our ap- 
proximation with reduced self-consistency (presented in 
Sect. is able to produce a propagating spin-wave-like 
mode, and the spectrum of auxiliary fermions shows a 
pseudo gap. In the fully self-consistent approximation 
(built in close analogy to FLEX), on the other hand, 
'spin- wave' mode and pseudo gap are absent. The related 
problem of strong collective Cooper-pair fluctuations in 
the normal state of 2D superconductors has been studied 
in the past^SiSLSS. The connection of these approaches to 
the present work is considered in Sect. Ivland Appendix 

El. 
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The order of the paper is as follows: The auxiliary- 
fermion formulation is introduced in detail in the next 
Section and the mean-field theory is discussed in 
Sect. mil . At zero temperature the mean-field theory re- 
sembles linear spin- wave theory'^^ (SWT). In Sect. IIVI 
the fully self-consistent diagrammatic approximation is 
presented. Static and dynamical properties are calcu- 
lated by analytical and numerical solution of the self- 
consistent integral equations. In Sect. the above- 
mentioned approximation with reduced self-consistency 
is explored in detail. It leads to results that qualitatively 
agree very well with what is known about the QHAF 
and the QNLtrM at low temperature. In Sect. IIVUI we 
apply the conserving-approximation method, and finally 
in Sect. IVll we estimate the fluctuations of the fermion 
charge Qi . The paper closes with a conclusion and some 
technical appendices. Calculations are performed at low 
temperatures T < 0.2 J , where scaling is observed. A 
possible crossover to quantum-critical behaviori^iiMSi at 
higher temperatures is not considered. Sections IIII Bi - 
ll V ClIVTl . and I VII contain material not directly related 
to the physics and can probably be skipped in a first 
reading. 



II. AUXILIARY-FERMION FORMULATION 

The constraint in Eq.(|3J) requires the fermion charge 
to be fixed to 1 at each lattice site i individually. This 
projection can be performed exactly in models like the 
Anderson impurity and Kondo model^is^SiiSi, which pos- 
sess just one atomic orbital with strong electron-electron 
interaction. In models containing a whole lattice of such 
orbitals the exact projection leads to a loss of the linked- 
cluster theorem, which prohibits infinite-order resumma- 
tion of the perturbation series^^i^ and self-consistent ap- 
proximations based on the skeleton-diagram expansion. 
The limit of infinite spatial dimension^S*^ has frequently 
been used to circumvent this problem. An alternative is 
to start from a mean-field theory**^, where the constraint 
is observed only in the thermal average. For the present 
case this reads, 



Q^^{Q^)={Ql)=J2(fUlc.)^l 



(6) 



and is introduced into the Hamiltonian jsj through a 
chemical potential /i-^ for the fermions. Due to particle- 
hole symmetry^ it is /i-'' = . The partition function 
can now be written as a standard coherent-state path 
integral^. 



with the action A = A° + A^- + A'^ , 

A° = dTj2J,^{T)[dr-flf]f^a{T) (8a) 

JO 

^ I, a 

A" = - Tdr Vh,(r)S,(T) (8b) 
Jo 

A' = j dTl^J,,S,(T)S,(r) (8c) 

with anticommuting Grassmann variables /, / , and [3 = 
l/ksT = 1/T . An arbitrary magnetic source field h has 
been added. The underlying lattice is a d-dimensional 
cubic lattice with lattice spacing a = 1 ; only the two- 
dimensional (2D) case d — 2 will be studied in detail. 

In the following the connected spin propagator is con- 
sidered, 



(9) 



^ (r.5f(r)5;(r'))-(5f(r))(5j'(r')) 

where fi = x,y, z denotes the spin component. It is cal- 
culated using a Feynman-diagram expansion^*^ with A"^ 
as perturbation. This will involve the Green's function 
of fermions, 

[G^Ar,r')^' - -{rrMr)fl,ir')) (10) 
= -(/.a(r)7,„,(r')) 
Time-r ordered expectation values are expressed via 



(M) = 



1 



Z[h] 



V[fJ] e-^M 



(11) 



Note that the fermion Green's function G is not an ob- 
servable quantity, since it depends on local U{1) gauge 
transformations^ of fia , //q . The susceptibility x j on 
the contrary, is built from the gauge-invariant physical 
observables Eq.jSj) . 

It is of great help that we can make use of standard 
Feynman-diagram techniques for a perturbation expan- 
sion in the interaction A"^ . This requires a mean-field 
like treatment of the constraint Eq.|(HI), since then per- 
turbation theory starts from (effectively) free fermions 
and Wick's theorem can be applied. In principle the 
treatment of the constraint can be improved by general- 
izing the chemical potential fi-^ to a fluctuating Lagrange 
multiplier*'^, but we restrict ourselves to the simplest 
approach^^. Nevertheless, in the approximations to be 
discussed below the pseudo fermions do not develop a fi- 
nite hopping amplitude connecting different lattice sites 

and anomalous expectation values of the form 



Z[h]^ / !?[/,/] e 



-A 



(7) 



ifiafja' 







are absent. As a consequence the Green's function of 
fermions is always local, i.e., 



G,. 



Sij Gi 



(12) 



Since unphysical expectation values are zero the local 
charge Qi introduced in Eq.JBJ is most likely conserved. 
Therefore spurious unphysical contributions to physical 
correlators like Eq.© are due to states with Qi ^ 1 
introduced by the trace in the thermal average Ea. (|ll|l . 
but these are not generated dynamically. We return to 
this issue later on in Sect. I VII . 



III. MEAN-FIELD THEORY USING 
AUXILIARY FERMIONS 

The Heisenberg term A"^ is a two-particle interaction. 
A mean-field decomposition of A'^ leads to an effective 
HamiltonianSii^SiSi^ with self-consistent order param- 

eters (AU""'/..') = 2(S.), , (Mfn) ■ 

These correspond to phases with different spontaneously 
broken symmetries, namely magnetic order, the so-called 
resonating valence bondS (RVB) or the flux phas o^^i^^ , 
and a d-wave paired state of the fermions^^iSiSij respec- 
tively. In the undoped case considered here, the latter 
two are equivalent by particle-hole symmetry'*^ . The flux 
and the paired state order parameters are unphysical, 
but may correspond to physical regimes observed exper- 
imentally (the pseudo-gap in the cuprates). The respec- 
tive transition temperatures are interpreted as crossovers. 
The magnetic order parameter J2a.a' ifL'^""' f^a') , on 
the other hand, is a physical (gauge invariant) quantity; 
a non-zero value indicates a true phase transition. If 
the free energies of these phases are compared in an un- 
restricted mean-field calculation, it turns ouli^i*^ that 
only (S) becomes finite in the Heisenberg model, indi- 
cating a Neel-ordered phase below some transition tem- 
perature Tn ~ J ■ This is consistent with the fact that 
the ground state of the Heisenberg model is close to Neel 
orderi. Only for finite hole filling in the i- J- model a 
phase is observed at low T , where AF and d-wave pair- 
ing order coexist22i. Therefore we do not consider flux or 
pairing instabilities. 

Note that the large-iV approachSSi^ leads to a dif- 
ferent ground state: The spin-index of the fermions is 
generalized from 2 to TV states, and the limit N oo 
leads to a saddle-point theory giving a flux phase with 
ifiafja) 7^ 0, while magnetism is suppressed. 



A. AF Order and Spin Waves 

With the source-field h set to zero, the Green's func- 
tion H10() . H12I) of the fermions reads 



Gi{iuj) — [iuj + fi-^ — Si(iaj)] 




n 





FIG. 1: Diagrammatic representation of E, D, A{S'^) in mean- 
field approximation, Sect. Illl'Xl . The thin dashed line is the 
bare interaction J in Ea. l|8c^ . the thick dashed line the renor- 
malized one Ea. l3UII . and the full lines are the Green's func- 
tion of pseudo fermions Ea. H18|l in mean-field approximation. 
Eg. 1151 . The dots represent Pauli matrices xl/2. A closed 
fermion loop requires a trace in spin space. 



uj = {2n + 1)tt/ P denotes a fermionic Matsubara fre- 
quency, G and S are matrices in spin space, and the 
chemical potential is always /i-^ = by particle-hole sym- 
metry. In mean-field the self energy Sj is given by the 
Hartree approximation depicted in Fig.n^tnd reads 

S,(z^) = ^J,, ^ia''iTrKG,(0,0+)] 

3 M 

Not shown in Fig.^are the (exchange) Fock and (anoma- 
lous) Bogoliubov diagrams, which correspond to the 
above-mentioned flux and pairing order parameters, re- 
spectively, and do not contribute. If we assume Neel 
ordering in z-direction on sublattices A and B, i.e., 
(5^) — —{Sg) , the self energy becomes 



'^^'"^ with rriA = ^ 2dJ (S-^) 



am A 



(14) 



The coordination number 2d has been inserted. The 
Green's function ((T^ now reads 



GA^iuj) 



tuj — (7 niA 



(15) 



and similar for B with iua — » ms = —niA ■ Inserting 
into (5^) delivers an equation for the Weiss-field m , 



~ 2 ^ 3 ^ iLO + am A 

a—±l iuj 

= ^tanh(!^) 
2 ^2T' 



(16) 



(13) 



Antiferromagnetic order 7^ appears below the 

transition temperature Tjq — jj ~ J/2. The aver- 
aged local density of states of the fermions, p{uj) = 



5 



W7 irlm^Tr [a'G,{Lo - iO+)] , becomes 



p{lu) = -[5{uj - itia) + S{lj + niA) 



(17) 



The spectrum has a gap 2mA as is known from the spin- 
density- wave state. However, our pseudo fcrmions do 
not show any dispersion, and hence p{ijj) is (5-hke. For 
tua = the bare spectrum p°(a;) = 5{ijj) is recovered. 

In the ordered state T < the magnetic excitations 
should be given by spin waves. In order to compare to 
the resuhs from hnear spin-wave theory, we calculate the 
spin propagator within the mean-field approximation. 
It has the RPA form 



([i + ./n]-in) 



(18) 



with the irreducible part 11 given by the fermion bubble 
shown in Fig. ^ 

nf^'(zi.) = -4 ^ ^Tiia^G.iiLo + iiy)a''"G,iiu)] (19) 

iuj 

V denotes a bosonic Matsubara frequency. Inserting 
Eq. (|15|l yields the transversal components 



2mA {S\) 
{2mAY - [wf 



(20a) 



and 



= -ny^{^^)^-^ '\f^]. (20b) 
[2mA) — yiv) 

—11^^ . The longitudinal 



Y\xx Y\^y 



response is exponentially small and can be ignored, 
~ exp(-/3|TOA|). The susceptibility ^ is 
transformed into wave-vector space using 



^'(q,q') = ^Y.^f 



'(q,q') = ^5]nr'e-('i-'i')^- 



(21) 



(22) 



The second line leads to 

n--(q,q') = n^^jq,q, , n-^(q,q') = n^Xq'±Q 

with the AF ordering vector Q = (tt, tt, . . .) and 6=*=*'^^' ~ 
+1, —1 for Hi ^ A, B . The interaction becomes 

J(q) = 2dJ7(q) with 7(q) = ^ ^ cos(g^) (23) 



Using Eqs.lEOI), lESl, with J(q±Q) = -J(q) gives 

1 - 7(q) 



X^^(q,*^^) = 2dJi{Sl)f 



f^2(q) _ [^^) 
iv 



(25a) 
(25b) 



f]2(q) _ (ivY 
These response functions have poles at ±ri(q) , 

r!(q) = 2dJ\{S\) I - 7'(q) - co|q - Ql (26) 

57 (q) resembles the magnon dispersion from linear spin- 
wave theorySS, (SWT). At low energy it has two branches 
with linear dispersion for q ~ and q ~ Q . For T <C TV 
the spin- wave velocity cq is temperature independent (up 
to exponentially small terms) and given by 



Co 



2\/d J I (S*!) I ~ 2\/d J S = V2 J 



This also reproduces exactly the result from SWT. The 
response function measured in neutron-scattering exper- 
iments is ^"(q,^) = Imx^^(q, + »0+) and comes out 
from Ea. (|25a|) as 



^ vi-7nq) 

X [(5(w-r2(q)) -(5(t^ + f7(q))] 



(27) 



The spin-wave excitations remain undamped (^-like) 
throughout the whole Brillouin zone, since max(r2(q)) = 
2dJ|(5^)| = 2|to^| at the magnetic zone boundary just 
reaches the gap 2|m^| for charge fluctuations of the aux- 
iliary fermions. The spectral weight of the magnons 
is measured by the transversal static structure factor 
5^(q,q') = - 5^(q)Jq,q, with 



5^(q) = - / dc.[l + .g(c.)]x"(q,c.) 
. meoth(-^(^^ l-^(^) 



(28) 



2r Vi-7^(q) 



g(Lo) denotes the Bose function. Eas. (|27|l . (|28|l result 
from mean-field theory and apply to any dimension d > 
1. At zero temperature KS*^)] coth(.)/2 1/4, and 
S^{q) becomes identical to SWT. 



B. 



1/5 Corrections 



By transforming Ea. (|18|l in q-space we find 

x""(q,q') = x^nq,q') = x""(q)'5q,q' (24a) 
x"''(q,q') = -x^"(q,q') = x"^(q)'5q,q'±Q (24b) 

and Ea. p8|) takes the form (omitting the q arguments) 



(1 + JW/) 

jn7 



(1 




The mean-field theory presented above can be obtained 
as the first approximation in an expansion in 1/5* . Fol- 
lowing Refs. I6QIi61m62[ the fermion is decorated with an 
additional 'orbital' quantum number k , ftcK , with 

K = 1, 2, . . . , iVs , A^s = S -M/2 for half-integer spin S . 
On each lattice site i there are now 2Ns — (25* -I- 1) states 
available. An expansion in 1/A^s ^ l/S* is performed by 
re-scaling the interaction J J — J/Ns with J = const. 
in each diagram and in Eq. (|18|l , and counting a prefactor 



Ns for each closed fermion loop. The limit Ns — > oo 
then produces the mean-field approximation. 

To complete the comparison of the auxiliary-fermion 
approach to spin-wave theory we also sketch the calcu- 
lation of the first correction ~ l/iVs to the staggered 
magnetization, (5f) = {Sf)^^^ + A{Sf) . Note that 
(5-2)^^ ^ {1/Ns)° . The l/iVs-correction is shown in 
Fig. n and reads 

^iS') = 23^EE^f''Mx (29) 

where we already returned to the physical case S — 
1/2 ^ Ns = 1, J = J ■ A renormalizcd spin interac- 
tion has been introduced, 



-Jij + E ■^'^ ^Ik i^'^) Jkj 



(30) 



l.k 



The local Di = Da appearing in Eq. (|29|l has no contribu- 
tion from the bare interaction, since Ju = . However, 
propagating spin excitations described by the suscepti- 
bility xik contribute to Di . With Eq. (|21|l it follows 



q.q' 



(q,q') 



(31) 



X j(q'; 



p«(q-q')R-i 



Ea. H29|l becomes for the A-sublattice, with the fermion 
propagator 1)15(1 inserted. 



4/32 E E X 

iuj,iu a— ±1 

aD'j^{iv) + iD''^{iv) 
[iLU + aniAYiii^ + iv — amA) 



(32) 



X^^ = X^^ s-nd x**^ = ^X^^ have been utilized. The 
longitudinal ^ 11^^ has vanishing weight and does 
not contribute. The transversal components read with 
Eq. (|31|l and the relations given in the text below Eq. (|21|l : 



J'(q)x""(q,*^^) 



= -^E'^'(q)^"'(q'*'^) 



T)XX 



r^xx 



'D^2' tlifi i?-sublattice. 



These 



are inserted into Ea. l|S^ . using the mean-field expres- 
sions H25f) for the susceptibilities. Furthermore Ea. H23|) 
and 7(q) — are taken into account, turning Ea. l(^ 
into an expression where the sums over fermionic (w) and 
bosonic {v) Matsubara frequencies still have to be per- 
formed. This is done with 

1 



^E^M 



i^f[z)T{z) 
J 27ri 

i^S(z^.) = -l^iicoth(/3z/2)6(z) 



n 




D 



+ 




+ 



+ 





+ 



FIG. 2: The self-consistent approximation for the paramag- 
netic phase at T > discussed in Sect. IIVI . Full lines de- 
note the fermion propagator G first introduced in Eg. 11311 . the 
thick dashed line is the renormalized interaction D defined in 
Ea. H^O^ . the thin dashed line represents the bare interaction 
J . Shown are the fermion self-energy E, Eq.lEHl, the irre- 
ducible bubble n , Eq. I|19|l , and the renormalized interaction 
D . Dots are Pauli matrices x 1/2 . In the paramagnetic phase 
the traces in spin space are easily performed, leading to the 
set of Eas. l|35|l . 



f{z) denotes the Fermi function, and the contour encir- 
cles the real axis, where !F{z),B{z) may be non-analytic. 
At low temperature (i \mA \ 3> 1 the result is 



A|(51)| = - 



25* 



(33) 



2N, 



E 



coth( 



a(q) 

2T 



Vi-7^(q) 



In the ground state at T = , where coth(il/2T) 1 , 
Eq.(|2SI) represents the reduction of the staggered mag- 
netization due to quantum fluctuations. It resembles ex- 
actly the l/S'-correction known from SWT. In 2 dimen- 
sions it is finite, with the numerical value A| (5*^)1 = 
-e = -0.197, i.e., for (S"^) > we have [(5^) + 
A(S'^)]/(5^) ~ 0.61. At T > 0, however, the diverging 
number of thermally excited magnons with low energies 
destroy magnetic order in 2 dimensions (theorem of Mer- 
min and Wagner): From Ea. (|33|l we get for wave vectors 
q close to or Q , where ri(q) < T and the thermal 
prefactor coth(f^/2r) ~ 2T/n , 



T > : A(5^ 



T 



reg. 



The integral diverges in the infrared k —^ for d < 3 , in- 
dicating the well-known breakdown of perturbation the- 
ory at any finite T in 2D. 



IV. SELF-CONSISTENT APPROXIMATION 
FOR FINITE TEMPERATURE IN 2D 

In order to fulfill Mermin- Wagner's theorem, i.e., to 
suppress the mean-field transition temperature down to 



7 




filling. Hence the fermion propagator Gi 



FIG. 3: Self-energy diagrams of ladder type for the exchange 
particle-hole ('RVB/flux') channel (top row) and the particle- 
particle (pairing) channel (bottom row). The contribution of 
jf^ and if has almost no effect on the results, see text, Sect. 

Hvl . 



zero, we seek an approximation where the susceptibility 
X is coupled back onto itself in a self-consistent fashion. 
Such an approximation has been proposed originally for 
the Hubbard model^'^, commonly referred to as FLEX. 
For the present model it takes the form shown in Fig. 
The irreducible part H in Ea. H18|l is approximated by a 
bubble of renormalizcd fermion Green's functions. The 
latter are in turn determined by x through the effective 
interaction D , which enters the fermion's self energy on 
one-loop level. When temperature is lowered and the 
mean-field transition at T ~ J is approached, the grow- 
ing susceptibility x should modify the spectrum of the 
auxiliary fermion and therefore H , such that the increase 
of X ^-nd accordingly the transition is suppressed. 

The irreducible part H and the fermion propagator G 
are given by Eas. lfT^ and l(T^ . respectively. The self 
energy depicted in Fig. |21 reads 



E(icij) 



Without external field and in the paramagnetic phase, 
h — (S) = , the fermions are spin degenerate, G ~ 
(j^G , S ~ a^Yj , and the response becomes isotropic, 



now simplify to 



-H. Equations (HHJ, (113, EOl, El 



Ii{iv) = - 
X{(\,w) = 



-y 



1 + ./(q)H(w) 



G(iw + w)G{iuj) (35a) 

(35b) 
(35c) 



E(jw) = ^^D{w)G{w ^ i^) (35d) 



Giiijj) = [iw — S(iw)] 



(35e) 



Note that the first bare J to D in Fig. [21 does not con- 
tribute, since Ja = or equivalently ^ J(q) = . 

It has been emphasized above that the non-local ex- 
pectation values ifJafja) and {fi^fji) are zero at 1/2 



Gi = G and the bubble H^- = Sij Hi H^ = H are local, 
i.e., wave- vector independent. The susceptibility x(q) , 
however, is q-dependent through the Heisenberg inter- 
action J(q) = 2dJ'-f{q) occurring in Ea. (|35b|l . The lat- 
ter has been introduced in Ea. (|23|l . However, the argu- 
ment for the fermions being local is based on mean-field 
theory. If fluctuations are included in the fermion self- 
energy, the ladder-diagram series shown in Fig. O should 
also be taken into account. While S in Fig. [21 captures 
the effect of magnetic fluctuations on the fermions, S 
and if^ contain fluctuations in the 'RVB/flux' channel 
and the pairing channel, respectively: A large numerical 
contribution from S would indicate a close-by insta- 
bility to a flux phase {{f} fj) ^ 0), whereas Tf^ would 
dominate if the pseudo fermions were instable to pair- 
ing {{fifj) 7^ 0). In the 1/2-filled case considered in 

this paper, both E are equivalent by particle-hole 
symmetry. The numerical calculation presented in Sect. 
lIVBI has been repeated with the two extra channels from 
Fig. [3 included. Even at the lowest temperature reached 
the results were only altered quantitatively; numbers re- 
ported in Sect. UVBl came out different by ~ 20% . The 
ladder sums (susceptibilities) appearing in jf^^^^ remain 
small, i.e., ~ 1/J . Apparently, the self energies from Fig. 
[31 can be omitted and will no longer be considered. 

In the following we calculate the antiferromagnetic cor- 
relation length ^(T) and the spectrum S'(q, cj) of mag- 
netic fluctuations from Eas. lp?^ . 



A. Static Approximation: Correlation Length, 
Energy Scale, Dynamical Scaling 

In order to proceed by analytical calculation we ap- 
proximate the effective interaction D in the self energy 
Ea. (|35d|l by its static value, D{iv) = D{Q)5ufi ■ This 
static (classical) approximation should be reasonable if 
the characteristic energy loq of spin fluctuations is small 
compared to temperature. This is the result of earlier 
studies mentioned in the introductory Sect. ^ and will 
also come out of the following calculation as well as the 
numerical solution presented in Sect. lIVBi b elow . 

With D{w) = £'(0)5^,0 Ea. l|35d|) turns into 

E(zc.) = l-D{Q)G{iLo) = (^) ' G(zc.) (36) 



where a coupling constant Uf — •\/3-D(0)//3 has been 
introduced. Ea. H35cl) requires the static susceptibil- 
ity x(qi 0) . At low temperature wave vectors q ~ 
Q = (tt, tt) close to the Neel ordering- vector Q give the 
strongest contribution, allowing to expand Ea. H23|) as 



J(q)~-2dJ+(q-Q)V 



leading to 

x(q,o) 



H(0) 



1 



1 



1 + J(q)n(0) Je-2 + (q-Q)2 



(37) 



(38) 
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where the antiferromagnetic correlation length has been 
identified as 



= 2d 



\2dJIl{0) 



1 



(39) 



The spatial dimension d has been written explicitly for 
clarity. In Ea. (|35c|l we may set J(q) ~ ^(Q) = ~2dJ , 
and it follows for the coupling constant 



fc2 



In 2 dimensions this becomes 
24 



(40) 



(41) 



where fee ~ 1 denotes a wave- vector cut off coming from 
the Brillouin-zone boundary. In higher dimensions c? > 3 
we get 



(42) 



The Dyson's equation (|35e() with Ea. H36|l inserted reads 
after analytical continuation iu ^ z 



1 



G(z) 



G{z) V 2 

At the real axis, z — u) + iO+ , this leads to 

G(. + zO,) = ^fc^ + 

2Q(\uj\-ujf) 



(43) 



w + sign(a;)^/ci;2 - (u;/)^ + iO+ 

In the atomic limit J = the spectrum of the fermions 
p(uS) = —^lmG{uj + iO-i-) is a delta peak at zero energy, 
p^{u!) = d{uj) from G'°(2) = Through the inter- 

action with spin fluctuations p is apparently broadened 
into a semi-elliptic spectrum with the coupling constant 
LOf appearing as a high-energy cut off: Ea. (|l^ gives 



piu;) = e{^f - 1^1) Jl-iiu/iOff (44) 

TTLUf » 

Eq. (|44|l is plotted in Fig. |S1 together with the numerical 
solution described in Sect. II V Bl . In order to close the set 
of self-consistency equations, 11(0) has to be calculated 
from Ea. (|35a|) . 

n(0) ^_ r de/(e)p(e)ReG(e) = —^T/ujf) (45) 

J-oc ^} 



with 



Ax 



27r 



37r 



The r.h.s. results from a Sommerfeld expansion of the 
integral. 

When temperature is lowered in the 2D system and 
the correlation length ^ grows, the fermion's cut-off en- 
ergy w/ increases (see Ea. (|41ll ). which in turn reduces the 
fermion spectrum p{Q) ^ l/iOf and the static irreducible 
part n(0) ^ . Therefore the antiferromagnetic sus- 
ceptibility x(Q,0) = n(0)/[l - 2djn(0)] is shifted away 
from the critical point, and the phase transition to AF 
order is suppressed self-consistently. In d > 3 the cor- 
relation length ^ does not enter loj , Ea. (|42|l . and this 
mechanism does not apply, allowing an AF transition to 
occur. 

An explicit expression for the correlation length is ob- 
tained from the solution of Eas. H39(l . H41(l . I|45(l . omitting 
terms of 0{^^'^) and 0{T'^) at low temperature. For 
d = 2 it follows from Ea. lP|l 



n(o) = i/u + oic"^) 



(46) 



This is inserted into Ea. (|45|) . which is then solved for ujf , 




UJf_ 

Finally, Eq.gTJ 
exp{(7r(..0)2)/(24rj) ^ 
order, 

e(r) cx exp(a, J IT) 



with LlPj 



gives the 
+ 0(r/J)}: 



with 



result 
and 

32 
27^ 



OTT 



in leading 



0.38 (48) 



Accordingly the correlation length diverges exactly at 
T = . The exponential behavior of ^ (T) is simply due to 

the density of states in 2D, M(yS) = J ^ 6{uj - J(q)) ~ 
0(4J— |(jj|)/ J , which is constant near the band edges and 
therefore leads to the logarithm in Ea. (|41|l . The corre- 
lation length (|48() is consistent with the theoretical work 
mentioned in Sect.Q]. The spin stiffness ps — asJ/2'K ~ 
0.06 J extracted from Ea. (|48|l comes out too small com- 
pared to the literature, where ps — 0.18 . The latter value 
is also the result of exact diagonalization studies in the 
AF ordered ground statai at T = . 

In addition to the correlation length, the dynamics 
of spin fluctuations is an important issue. The damp- 
ing of spin fluctuations is induced by the imaginary part 
of n(a;) , i.e., the spectrum of particle-hole excitations. 
From Eq. (|35a|) we get at T , with the smooth fermion 



spectrum (|44f) 

Imn(tj + 2:0+) = 



de p{lo — e)p{—e) 



'-piQfu 



The real part is approximately calculated using 

z0+) 



Ren(w) 



1 . ^Imri(e, 

TT J e-L 
n 

de- 

-n £■ 



OJ 
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where a high-energy cut off fl J for spin-fluctuations 
has been introduced. f2 can be fixed by the condition 
Ren(O) = n(0) , and it foUows 

Ren(w) = n(0) [l - [oj/nf] + 0{uj'^) 

with n = n(0)/p(0)2 = fJ. Here n(0) = 1/4 J and 
p(0) = l/7rw° = 3/16J from Eqs.^, ^ have been 
used. The susceptibiUty Ea. (j35b|l at the real axis iiy 
LJ + iO+ becomes 



n(w + iO+) 



x(q) n(o) 



^(q) 



'2 (ri) 



f- 

4 \n 



with the static susceptibiUty x(q) stated in Ea. (|38|l . In- 
troducing a spin-wave velocity c and -damping constant 



7, 



4n(o) jn^ 



(1.47 J)2 , 7 



1.90 J 



the dynamical susceptibility takes the familiar form 



r(q) ^ Jk^ . This is in strong contradiction to the above- 
mentioned theories, where for k ^ the structure fac- 
tor features two peaks with some dispersion w — ±f2(q) 
and finite width. These peaks are interpreted as spin- 
wave like excitations, visible only at short enough wave- 
lengths <C ^ . 

Since the w^-term in Ea. l49|) does not contribute, 
the dynamical scaling property (jS)), Q is almost triv- 
ially fulfilled: From Eqs.lgHI), (EH) it follows S'(q,Lj) = 
^^Im [1 + (fc^)^ — t(cj/(jJo)] 7 with the energy scale 



c 

^0 = r(Q) = — (X JC 
7r 



(52) 



The static structure factor results from 5**(q) 
i/dc.5(q,c.) rx(q,0) = ^[1 + (fc^^j- 
ducing Eqs.©, with the scaling functions 



repro- 



ip{x) 



1 



$(x,2/) 



(1 + x^Y -I- 



(53) 



For any fixed wave vector x = k^ the 'shape function' 
$(a;, y) , y = w/wq is a simple Lorentzian located at w = 
, normalized to $(0, 0) = 1 . The inset of Fig. 01 shows 
$ from Eq. H53I) , the main figure the corresponding result 
from the numerical solution. 



X(q,C.) = y [c2(fc2+^-2)„^2_^^^] (49) 

where k = (q — Q) . 

The calculation seems to yield spin-wave like excita- 
tions for wave vectors A: ^ far enough off the Neel 
vector Q , with the familiar dispersion f2(q) = ck . How- 
ever, the damping 7 appears much larger than the spin- 
wave energy, i.e., the 'spin waves' are totally overdamped: 
For wave vectors well inside the 1. B.Z., k ^ 1 , it is 
\uj\ ~ ri(q) ^ Jk <^ J and thus 7 ^ |a;| . Therefore the 
Lo^ in Ea. (|49|l does not contribute, and 

X q, w = X q^T^rr^ , F q = — - — 

The dynamical structure factor 

Imx(q,w) T 



(50) 



S{q,u;) 



1 



imx(q,w) 



(51) 



shows a single relaxation-like peak at a; = . At long 
wavelengths A: <C the half-linewidth of this peak is 
F(q) cx r(Q) J£,~^ , which turns to zero as the AF 
transition is reached (critical slowing down), with the 
dynamical exponent z — 2 . This is compatible with 
hydrodynamic theory^- and the studies of 5'(q, w) ref- 
erenced in Sect. []. Our result z = 2 , however, differs 
from the z = 1 found in these studies. At short wave- 
lengths, k ^ , the structure factor still consists of 
a single peak at a; = , with a fc-dependent linewidth 



B. Numerical Results 

So far we have considered a static approximation, 
where the susceptibility Ea. lj35b(l has been taken in the 
classical limit, x(qi*f^) = x(qiO)'^i',o inside the effective 
local interaction D , Eq. (|35c|) . The set of equations H35|) 
has also been solved numerically, taking all frequency 
dependencies into account. Eas. H35|l are re- written at 
the real axis, see Ea. ljAl|) in Appendix lA II . and solved 
by numerical iteration. A stable solution is obtained 
for temperatures 0.5 J > T > 0.054 J, corresponding 
to correlation lengths 0.7 < ^ < 101 , in units of the 
lattice spacing. ^ is extracted from the static suscepti- 
bility at q = Q using Ea. H39|) . it is well reproduced by 
^(T) = const. X exp(a J/T) with a = 0.30 . The fit of the 
numerical data ^(T) is performed by a linear regression 
of [iln(^)] to the function [a + t\n{const.)] , t = T/J. 
The exponent = 0.38 from the static approximation, 
Ea. (|48() . almost agrees with the numerical result. 

The numerical solution is also consistent with the dy- 
namical scaling hypothesis: If the static scaling function 
(p{x) — S"**(q)/5"'*(Q) is plotted against x = kS^ , curves 
lie on top of each other for temperatures T < O.IJ where 
f > 10 . (p{x) is well described by the simple Lorentzian 
shape ()53|l of the static approximation. The energy scale 
ujQ is extracted from the numerical output via 



g^*(q = Q) 
S{q = Q,Lu^O) 



(54) 
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4 6 
Energy co/co^j 

FIG. 4: Main figure: Dynamical structure factor from 
the numerical solution of the self-consistent approximation 
Eas. (l35t . Plotted is the scaling function <^(k^,uj/u>o) at 
low energies ~ ljq • "l" is shown for fixed wave vectors 
k — 0,^~^, 2^~^, 5^~^ , each for 2 different temperatures cor- 
responding to the correlation lengths ^ = 12.1, 101.2 . The 
curves lie almost on top of each other, indicating that dy- 
namical scaling holds for ^ > 10 . Inset: Analytically calcu- 
lated scaling function <&(fc^, cj/tJo) , Eg. 1531 . from the static 
approximation. 



which gives wq up to a constant prefactor, which is chosen 
such that $(0,0) = 1. The scaling function (^{x,y) is 
then gained from 



(55) 



and plotted against x = , y = lu /loq ■ Fig. ^ shows 
^{x,y) for two temperatures, where C ^ 10 • Apparently 
scaling is well obeyed, and the shape of $ agrees with 
the static approximation. In particular, it consists of 
a single peak even for k£, > 1 . In Sect. a modified 
self-consistent approximation will be presented, where 
^{x,y) features spin-wave like propagating excitations 
for fc^ > 1 ; also the dynamical exponent will come out 
correctly as z = 1 . 

For completeness the spectrum p{uj) of auxiliary 
fermions is shown in Fig. [S]. The delta-peak of the bare 
{J — 0) fermion is broadened into a continuum of width 
^ 2J, qualitatively similar to the static approximation. 
Eg . (|44|l . Further comments on p{u!) will be made in Sect. 
CVDl. 



C. Conserving- Approximation Approach: 
Vertex Corrections 

We now turn to the so-called conserving- 
approximation method, which has been left aside 
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FIG. 5: Spectral function p{io) of the fermion propagator 
Ea. H35e^ . Full line: numerical solution of Eas. ()85^ . Dashed 
line: static approximation, Ea. l|44^ . The arrows indicate 
delta peaks, representing p{uj) in the antiferromagnetic phase 
at T = in mean- field theory, Eq.lUZJ . 



so far: The spin susceptibility (j^ is given by 



(r') 



-S^W[h] 



h=0 



5/if(r)5< (r') 



h=0 



/i, /i' — x,y, z . h'^ are the components of the mag- 
netic source field h . is the generating functional 
(free energy) W[T,h] = - In/ X>[/,7] c"^ to be cal- 
culated from the action ||SJ) of the model. W can be 
expressed as a functiona l^'^i^^ of G and E , M^[r, h] = 
$[G]-2Tr[SG']-2Trln[-Go ^-HE] . The functional $[G] 
appearing here is related to the fermion self-energy via 
E[G] = i(5$/(5G. The existence of a ^-functional for 
a particular approximation of S is a sufficient condition 
for X to respect conservation law a^^i^^ . The self-energy in 
Fig. 13 Ea. H35|) can be derived from a ^-functional similar 
to the FLEXS. 

X is given by the same RPA-formula H18|l as x , but 
the irreducible part 11 from Fig. |21 has to be replaced 
by n shown in Fig. If the theory was exact, 11 and 
n would be identical. In any approximation, however, 
n contains additional vertex corrections, and x differs 
from the correlation function x of self-consistency 
equations, e.g., Eas. (|35|) . If now x diverges, indicating 
a magnetic instability at some finite temperature Tjv , 
this is not coupled back into the self-consistency cycle. 
As a consequence the spurious phase transition is not 
suppressed, and x becomes invalid below Tjv . 

To elucidate further on this important point, we cal- 
culate the correlation length ^(T) from the susceptibility 
x(q, 0) . We employ the static approximation from Sect. 
IIV Al using the fermion propagator H43|) and the self- 
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FIG. 6: Top: Irreducible bubble 11 including the vertex corrections from the conserving approximation-scheme. Bottom: 
Bethe-Salpeter equation for the vertex function A . Lines are those introduced in Fig. |5| calculated self-consistently from the 
corresponding equations l|85|l . 



consistent D , Eq.JSSfijl, as input. The bubble including 
the vertex corrections is shovi^n in Fig. Eland reads 



G and A denote the renormalized fermion propagator and 
the vertex function, respectively; both are matrices in 
spin space. Whereas the former is always local, the latter 
in general depends on lattice sites i,j . In the disordered 

phase (S) = we have Gi = Gcr° and ftf/' = Hf/p^' . 

For the correlation length we set w = , and IT reads in 
wave- vector space 

n^(q,0) = -i^GM^A^qjzc.) (56) 



A''(q;*c.) 



lTrKA^(q;zc^)] 



The Bethe-Salpeter equation for A is depicted in Fig. El 
for ii/ = it takes the form 

F(q;*c.) = ia^+ 



1 1 

4^ 



D{0)G{iujf ^ a'^A''(q;iw)a''- 



fi=x,y,z 



^Y.G{iu;)'TT[a^A^{q;iCoy^'] 



D(k,ii^) denotes the renormalized interaction Ea. H30() 
and D{iv) its local version Ea. H35c|l . In the static ap- 
proximation set D{iv) — D[Q)5u,o ■ The last term in the 
equation above corresponds to the last two diagrams in 
Fig. 0. In the static case ii^ — considered here, these 
diagrams cancel, since (a'^a'^ + ct'^ cr'^) — 2o'°(5^^' , and 
therefore their contribution to A*^ is (x Tr[(7^cr''] = . In 
addition, in the undoped case of the Heisenberg model it 
is G{iojy — , I = 3, 5, 7, ... , due to particle-hole 
symmetry. Thus only the first two terms in A contribute. 



giving a local vertex function 
1 1 



A^(q;zc.) 



4/3 



D{0)GiiLuy X^'iq-^iuj) (57) 



We used Tr[cr'^cr^cr'^A ] = (±) Tr[cr^A^] with (+) for ft = 
/i and {—) ior jl ^ fi . From Eas. (|57|) and (|56|) it follows 

n^(q,o) = n(o), 



n(o) 



2tt 



de /(e) Im 



G(e + i0+)-2 + i(u;//2)2 



with the coupling constant ujf — \J2>D{Q) / (3 introduced 
below Ea. H36l) . For the fermion propagator G the self- 
consistent solution Ea. (|43|l has to be inserted, resulting 
in 



n(o) 

m 



1 

^/ 

3 

TT 



(58) 



dx 



r:\/r 



c2 + 1/3 e^/* + 1 



37r 



and $(0) ~ 0.497. Taking the T-dependent cut off Uf 
from Eq. (|47|l , the vertex-renormalized bubble becomes 



n(o) 



i(0)-(r/a;?)'5(3-^$(0)) +0{T^) 



= - [(0.293) - (0.462) [T/jf 



Now the static susceptibility x is formed similar to 
Eq. (|38() , and the corresponding correlation length is given 
by ^ = Jx(Q,0) = jn(0) / (1 - 2djn(0)) , leading to 



T, 



1/2 



N 



T -T, 



with Tat = (0.305) J (59) 



N . 



The spin susceptibility x derived from the conserving- 
approximation principle contains a spurious transition to 
AF magnetic order at some finite temperature ~ 
0.3 J . The mean-field transition T^v = 0.5 J is only 
slightly suppressed. Close to Tjv the correlation length 
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^ shows a mean-field like behavior, rather than the ex- 
ponential one of the 2D system. The vertex correc- 
tions included in 11 apparently have a dramatic effect: 
n(0) , the irreducible bubble entering Eas. (|35|l . is forced 
by self-consistency to be constant, J(Q)n(0) = (— 1) -I- 
0{e~'^/^) , see Ea . (14611 . Accordingly the fermion cut- 
off cj/ , Ea. (|47|l . comes out such that this 11(0) is repro- 
duced if w/ is inserted into Ea. (|45|) . in particular the T- 
dependence oiojf cancels the one from <i>(i) , up to terms 
©(e"'^/'^). However, if the same is inserted into 11(0) , 
Ea. (|58|l . this cancellation no longer takes place, leading 
to the qualitatively different result for ^ compared to f . 

In conclusion, the susceptibility derived from the prin- 
ciple of conserving approximation leads to very unsatis- 
factory results. It is expected that the striking difference 
between ^ and f is not overcome if a more sophisticated 
approximation for the self-consistent 11(0) is used, since 
again 11 will differ from 11 . Therefore we will drop the 
conserving-approximation scheme in the following and 
take the self-consistent correlation function x &lso as an 
approximation for the spin susceptibility, x X • Since 
the sufficient condition for x to respect conservation laws 
is now lost, this issue has to be discussed separately. 



D. Absence of Short- Range Order: Overdamping 

A self-consistent approximation, designed in close anal- 
ogy to FLEX, provides the required suppression of the 
mean-field transition temperature Tjv — > in two dimen- 
sions. The correlation length ^(T) comes out qualita- 
tively correct. However, the dynamical structure factor 
S'(q, <jj) does not show any damped spin- wave peaks ex- 
pected far off the Neel vector Q = (tt, tt) . 

The lack of peak dispersion in S can be traced back 
to the fermion density of states p{uj) , displayed in Fig. 

There is no sign of a suppression around lo — Q even 
as T ^ and therefore any magnons are completely 
damped out. In the AF ordered state at T = the 
fermion's spectrum must show some sort of (pseudo) gap, 
allowing for well defined spin waves. In mean-field theory, 
there is even a true gap, see Ea. (|17|) . In the following sec- 
tion a modified (actually simplified) approximation will 
be discussed, where a pseudo-gap like suppression around 
w = is observed in p{uj) . In that approximation the 
spin dynamics come out quite satisfactorily. 



V. MINIMAL SELF-CONSISTENT 
APPROXIMATION (MSCA) FOR T > 

The most important property is the self-consistent 
treatment of the susceptibility x • This was achieved 
by renormalizing the propagator G of the fermion self- 
consistently, with x appearing in the self energy to one- 
loop order, see Fig.|21. In that scheme, the fermion propa- 
gator is subject to a self-consistency loop even if x is held 
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FIG. 7: The minimal self-consistent approximation (MSCA) 
for the disordered phase at T > discussed in Sect. IVl 
Eas. 160^ . The full line denotes the fermion propagator 
Ea. l60eL which contains the self energy E, Ea. ll60d|l . The 
thick dashed line is the renormalized interaction D, Eg. lGUcH . 
Note that in contrast to the fully self-consistent approxima- 
tion shown in Fig.|21 the fermion propagator in the self energy 
E as well as one of the fermion lines in the irreducible bubble 
n, Ea. l60aL remain unrenormalized (bare). These are given 
by G'^{iuj) = \/iuj . In the MSCA the only self-consistent 
quantity is the susceptibility x , Ea. (l60b| . that enters D . 



fixed, which smears out any pole structures. Fig.Qshows 
an approximation, where the amount of self-consistency 
is minimized, which we refer to as the 'minimal self- 
consistent approximation' (MSCA). In Fig. [7| only the 
susceptibility (the bubble series entering D) is a self- 
consistent quantity. The auxiliary-fermion propagator, 
which is itself not a physical object, is kept unrenormal- 
ized. The H shown in Fig.[71is a bubble of bare fermions, 
with D ~ X inserted to infinite order in the simplest way 
that preserves the analytic properties of H(z) . 

The set of equations corresponding to Fig. [7| actually 
is a simplification of Eas. (|35|) . With the bare fermion 
propagator G^{iuo) ^ l/iuj they read 



D{iv) 
G{iuj) 



VG'(^w-^w) — 
2 (j ^-^ ni.) 



1 + ./(q)n(ii/) 



4 B — ' lu) + IV 

iv 



(60a) 

(60b) 

(60c) 

(60d) 
(60e) 



Again we first consider the static approximation, i.e., let 
D{iv) = i:>(0)(5^,o in the self energy Ea. (|60H|l . This leads 
to Yj(iuj) — (clis)^ / and thus 



G(iu) 



1 
2 

2 



1 



1 



(61) 



4/3 



I?(0) 



That is, the fermion spectrum consists of two (5-peaks at 
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±u)s ~ J , which is much more close to the anticipated 
(pseudo) gap than the continuous one Eq.|05 . Within 
the static approximation the gap comes out as a hard gap 
of 2[Ds , it almost resembles the mean-field result for the 
ordered state at T = given by Eq. iQT)) . From Eq. (|60af) 
we get 

Here also the similarity to the mean-field expression 
Ea. (|20al) is apparent. At the real axis ii^ ^ uj + iO+ , 
we find using J(q) ~ —2dJ + Jk^ , k = q — Q , 

X(q, u;)^j [c^iC^ + e) -u?~ ic^O+] (62) 

with (? = Jw^n(O) and from Eq.ljSHl) . Omitting 
terms of ©(e-'^^O and Oi^"'^) , for d = 2 it is n(0) = 
1/4 J from Eg. ipi) . n(0) ^ l/4a;s , and thus cj^ = J. The 
latter has to be equal to Cos — w//2, with ojf already 
derived in Ea. (|41|) . The solution of the resulting self- 
consistency equation gives a correlation length and spin- 
wave velocity 

i{T) cc exp(a, J/T) , a, = J ~ 0.52 , c = ^ (63) 

6 2 

Accordingly, in the static approximation, the MSCA 
delivers sharp 'spin-waves'. Their dispersion ri(q) — 
y/ + c^k^ > A features a gap A — c^~^ , since at 
T > the spin-rotation symmetry is not broken. In 
the static approximation the dynamical scaling property 
©, Q is almost trivially fulfilled by Ea. (|^ . since the 
linear term ~ iuj does not contribute: From x(q, uj) = 
j[l + {k^y - {uj / A)'^ - iujO+]-'^ the structure factor ifH)) 
immediately takes the form ©, with the scales ^ and 
luq ^ A . The dynamical exponent comes out correctly 
as z = 1 . The infinite lifetime of spin waves originates 
from the hard gap 2ws = 2J in the fermion propagator 
(|61|l . If the static approximation in E is relaxed, this 
gap will be partly filled, and spin excitations acquire a 
finite damping. In that case both terms ~ iuj and ~ 
contribute in x(q, w)~^ , and dynamical scaling becomes 
non-trivial. 




Temperature T/J 



FIG. 8: Main figure: Correlation length ^(r) of the mini- 
mal self-consistent approximation (MSCA). Shown is the data 
from the numerical solution of Eas. l)6()^ and the function 
^(T) = 0.362^ exp(0.359 J/T) resulting from a linear regres- 
sion of the data, plotted as [t(ln(^) + ln(f))] , to the function 
[a + t\n{b)],t = T/J for T < 0.18 J. Inset: Energy scale 
u;o{T) of the MSCA, extracted from the numerical solution 
of Egs. lfciUl using Ea. (l54ll . Shown is the numerical data for 
ciJo(r)^(T) and the function tJoC = 0.690 T, obtained from a 
linear regression for T < 0.18 J . The extrapolation to T = 
is (.Jo (0)5(0) = 0.00003 J. 



[3 it is the particle-hole channel (magnetic susceptibil- 
ity). In both cases the ladder element (or bubble) is 
formed with one bare and one renormalized Green's func- 
tion. Apparently, in order to obtain a pseudo-gap in the 
fermion spectrum due to strong collective fluctuations, 
an approximation with 'reduced self consiste ncy ' is de- 
sirable. This has also been observed in Ref. |3a for the 
Hubbard model in 2D close to the antiferromagnetic tran- 
sition. 



A. Relation to Earlier Work 

Our minimal self-consistent approximation bears sim- 
ilarity t o a schem e that has been introduced early on 
in Refs. I36l3i67l and later used in Refs. I38l68l to study 
pseudo-gap effects caused by strong (resonant) Cooper- 
pair-amplitude fluctuations in low-dimensional systems. 
The physics is similar, except for the type of long-range 
fluctuations responsible for the formation of a pseudo 
gap. The fermion self-energy used in, e.g., Ref. |3^ con- 
sists of particle-particle ladders (pair susceptibility) con- 
voluted with a bare fermion propagator, whereas in Fig. 



Supposedly, the MSCA cannot be obtained from the 
stationary point of an approximate free-energy functional 
(<i>-derivable approximatioriS^iS&) , since both bare and 
renormalized Green's functions enter the self energy^. 
However, exact equations of motion in combination with 
the use of cumulants (connected Green's functions) pro- 
vide an alternative systematic approach. The basic ap- 
proximation is to ignore the effective interaction (ver- 
tex function) of 3 and more fermions. That technique 
has been introduced in Ref. |63 and applied to supercon- 
ductors in Refs. I3I37I . In Appendix ini we re-derive the 
MSCA for the Heisenberg model from the equations of 
motion. 
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FIG. 9: Main figure: Dynamical structure factor from 
the numerical solution of the MSCA, Egs. 11001 . at low en- 
ergies ~ . Plotted is the scaling function oj/tJo) 
as defined in Eg. 15511 . $ is shown for fixed wave vec- 
tors k = 0, ^"'^j 2^"^, 5^"^ , each for 3 different tempera- 
tures corresponding to ^ = 3.6, 44.4, 403.4 . The curves for 
^ = 44.4, 403.4 are almost indistinguishable. For wave vec- 
tors fc > a propagating 'spin- wave' mode becomes visi- 
ble. Inset: The scaling function $(fc^, tj/tJo) of the analytical 
calculation described in Sect. IV CI Eg. 1851 . for p = 0.25 . 
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FIG. 10: Main figure: Numerical spectral function p{uj) of 
the fermion (full line) in the MSCA, Eqs.||Snj, in the high- 
energy range ~ J . The dashed line is the high-energy ap- 
proximation Eg. 1771 from the analytical calculation. Inset: 
Numerical p{u-i) at low energies ~ tJo , shown is the scaling 
function p{ijj/ujo) = ^^p{^) introduced in Eg.ffH . The 
curves for ^ = 44.4 and 403.4 (dotted and full line, respec- 
tively) are almost identical. The dashed line is p{uj/ujo) from 
the analytical calculation. Eg . 1741 . 



B. Numerical Results 

In the following we solve Eas. HSUI) numerically, tak- 
ing the frequency dependence of D{iu) into account. 
Eas. H60|l are analytically continued to the real axis, iv — > 
u) + iO+ , leading to Eas. (|A2|l . These are solved by nu- 
merical iteration. Numerically stable solutions have been 
obtained for temperatures T > 0.08 J, corresponding to 
correlation lengths ^ < 400 lattice spacings. Correlation 
length (,{T) and energy scale oJo{T) are extracted from 
the numerical data using Ea. H39|) and Eg. 1)54(1 . respec- 
tively. For T < 0.2 J , where ^ > 10 , these are well fitted 
by 

C(T) cc '^eMaJ/T) , a = 0.36 , Uo{T) cx j (64) 

Fig. |S1 shows the numerical data and the fits. 

The result for ^ contains a factor 1/T in front of the 
usual exponential. In that respect ^ is different from the 
static approximation discussed above. The 1/T prefac- 
tor has been obtained for the QNLcrM , as long as the 
RG calculation is restricted to 1-loop ordcr^°. The en- 
ergy scale cjQ ~ I which we confirm by analytical 
calculation later on, is now consistent with a dynamical 
exponent z = 1 . The MSCA also features short-range or- 
der in the dynamical structure factor S'(q, uj) : The scal- 
ing function $(fc^,Ci;/wo) , extracted from the numerical 
data using Ea. (|55|l . is plotted in Fig.|^. Apparently the 



single peak at = splits into two peaks located sym- 
metrically to a; = , if the inverse wave vector is smaller 
than the correlation length, fc^ > 1 . (Negative energies 
are not shown in the figure). These peaks move away 
from cj = with increasing k , similar to damped spin 
waves. According to the analytical calculation, their dis- 
persion w(fc^) becomes Hnear for fc^ 1 , |a)|/a;o = 2fc^ . 
The numerical data obeys scaling for temperatures below 
^ O.IJ, where ^ > 10 . Dynamical scaling is also visible 
in the spectrum p{uj) of the fermions. The inset of Fig. 1101 
shows the scaling function p{w/u;o) = ^^p{uj) cx S,p{uj) 
at low energies ~ luq ■ The curves for ^ ~ 40 (dotted line) 
and for ^ ~ 400 (full line) lie on top of each other. The 
definition of p is a result of the analytical calculation, 
Eq. (|74|l . The main Figure ^| shows an overall view of 
p{llj) , the anticipated pseudo-gap structure is apparent, 
with two broad peaks near ± J . 



C. Analytical Calculation: Origin of Scaling and 
the Propagating Mode 

In order to confirm the numerical result for the en- 
ergy scale, luq (X T/^ , and to obtain some insight into 
the origin of the dynamical scaling numerically observed, 
we finally discuss an approximate analytical solution of 
Eas. H60|l for T <^ J and ^ ^ 1 , which does not make use 
of the static approximation. 

We start from the following set of equations, which is 
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derived in Appendix from Eas. H60() . 

e 1 



J 1 + (CA:)2 - R{uj) 
R{lo) = 16JC^n(w) 
n(w) = i^tanh(cj/2T)p(tj) + 

Z-"", ^ tanh(e/2r) 



= Ini[a; — 

o poo 

= — / d£lm£i(e) 
D{u) = — I dx 



coth(e/2T) 
w — e + iO+ 

1 



1 + x- R{lu) 



(65a) 

(65b) 
(65c) 



(65d) 
(65e) 

(65f) 



The correlation length ^ is given by Ea. H39|) as before, fc^ 
denotes a wave- vector cut off of 0(1) . The static ui = 
value of n required for Eq. reads 



n(0) - ^ / de tanh(e/2r)^ 
2 Jo £ 



(65g 



In Ea. (|65af) the static properties of x have been sepa- 
rated from its dynamics, which are determined by 11(0) 
(through £_) and Il{oj) (through R{uj)), respectively. See 
Appendix El for details. 

For the dynamical part R{uj) we consider the ansatz 



R{uj) = iiv/ujo +p{ijj/uJoY 



(66) 



ujq and p are independent of lo . This form of R{uj) , 
inserted in Ea. H65a|l . is able to describe the relaxation 
regime A;^ < 1 as well as propagating 'spin-wave' modes 
for fc^ > 1 . 

With Ea. (|66(l the imaginary part of the renormalized 
interaction (|65f|) becomes 

4J 

lmD{LL>) = — arctan(a;/cJo) ©(^c ^ (67) 

TT 

Iml?(a;) is essentially the local spectral function for spin 
excitations. In Ea. (|67|l all features of Iml?(a;) for u; > uiq 
have been replaced by a constant; loc ~ 2J is a high- 
energy cut-off. 

The energy scale loq is exponentially small, loq ~ J/S, ■ 
This already follows from the real part at w = , which 
results from Eq.jeSJ) as ReL'(O) = ^ M^cO and should 
roughly compare to the result of a Kramers-Kronig trans- 
formation applied to ReD{0) ~ ^ ln(cjc/^o) ■ 

To calculate ujq and p with Eq. as input, Eq. I|65c|l is 
expanded up to 0{uj'^) and the parameters are identified 
as 



ujQ T 



de 



tanh(e/2r) p{e) 



(68) 
(69) 



The energy scale ujq is extracted easily with the ob- 
servation that the self energy at w = is imaginary, 
Rel](0) = . This comes from the particle-hole symme- 
try of the Hamiltonian JSJ) at 1/2-filling. With Eas. l|nSe|) 
and l|67() we get lml](0) = —'SJT/ttujo , leading to 
p(0) = -(7rImE(0))-i = ujq/SJT . With Eq.JHHl a self- 
consistency condition for loq follows, with the solution 



Wo 



Apparently the dynamical exponent is z = 1 . 

The parameter p requires an integration over the 
fermion spectrum p{e) . Therefore we approximate S(w) 
for low w < cjQ <C T and high energies uj loq separately, 
at low T ^ J, ojc ■ For uj < ujq we have from Eq. (|65e() 
and lEZl, 



ImE(w);o,u 



3JT arctan(ti;/cL;o) 



(71) 



The real part is obtained from the Kramers-Kronig 
transform 



1 

ReE(t^) = - / de 

^ J-oo 



(72) 



as 



Rel](w) 



6JT 



low 



de 



arctan(e/a;o) 



3JT ln[l -I- (uj/LUoy 



2tt 



(73) 



Since lu ^ too the self-energy is large, |S(a;)|;oiu ~ 
JT/ujq oc 3> w , and the spectrum Ij65dp is given 
by piuJ)lo^u = ilm[l/E(cj)] + ol^-^) . Accordingly the 
fermion spectrum at low energy takes the form 



low 



3JT 

UJQ 

3JT 



Im- 



y 



i ln(l + 2/2) — i arctan(?/) 
Piy) : y^uj/ujo (74) 



where some scaling function p{y) has been identified. For 
?/ — > the previous result ifTHjl for p{Q) is recovered. p{lo) 
shows an incomplete gap, the spectral weight in the re- 
gion Lo loq IS 1/ JS^ . The scaling function p{uj/ujq) 
is plotted in the inset of Fig. ^| it qualitatively agrees 
with the result from the numerical solution. 

In the regime of high energy uj luq the imaginary 
part of S reads 

3 / 

lmJ:iio)h^ = -—cot^\iu\/2T)e{u;c-H) (75) 

and the real part is calculated using Ea. H72|l . 

coth(£/2r) 



Re2:(u;)ht = —Lu / de ■ 



3J 

L 

2tt 



de 



2T 



de- 
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A lower integral boundary ujq has been introduced, since 
in Ea. (|65el) the arctan from Ea. (|67l) cuts off the 1/e di- 
vergence coming from the coth(£/2r) . In the result- 
ing HeEhi we keep only the terms ^ l/uj showing the 
strongest u; dependence, that is, 



ReE(w) 



hi 



Ul 
LU 



3JT 



HT/lo^) (76) 



Thus we arrive at the fermion spectrum in the high- 
energy regime. 

This spectrum has two peaks ai uj — iiD , similar to 
the result from the static approximation appearing below 
Eq.lEOl, and LOa are of the same order of magnitude 
~ J . However, here the peaks are no longer (5-like but 
show a significant width. p(uj)hi is plotted and compared 
to the numerical solution in Fig. 1101 . 

With Eas. H74|l and H77|l the parameter p can now be 
estimated: The integral in Ea. H69|) is split into three re- 
gions and the tanh(e/2r) approximated. 



12 



-JT-' 



1 

2T 



~2 2 



2T 



de 



+ 

P{£)hi 



12 o , 

IT 



It + Ij] 



(78) 



The result 170(1 for has been inserted, a is a num- 
ber, awp indicates where piow and phi roughly match. 
We assume that a can be chosen constant at small T , 
where luq ^ awo ^ T ^ J, wc ■ We seek the leading 
T-dependence of p . 

With Eq.(|73|l the low -energy contribution becomes 



Iu>o = 1™ lo dxp{x)/{x'^ - y2) , that is. 



(79) 



For estimating It the spectrum 1(77(1 is expanded for < 
T <C J, (I" , i.e.. 



p{uj) 



hi 



— Im 

TT 

3JT 

and we arrive at 
3J 



ZUJ UJ 



(80) 



It = 



ln( — ) -ln(a) 

UJQ 



l + 0(§) 



The leading T-dependence of w is a constant, u) ^ J (see 
below), and we find 

/t cxl/ J^T (81) 



In the high-energy contribution Ij we use Ea. ((77() with 
coth(|cL'|/2T) ~ 1 . The dominant contribution to the 
integral in /,/ comes from the lower bound T . That is. 
Phi may be expanded as above, omitting terms 0{lu/J) , 
which results in p{uj)hi = (3 J/4[D'*) . With lo = 
const. ^ J its contribution to p is 



I J (X HJ/T)/J' 



(82) 



Collecting Eqs.ljTU, (|HU, together, the leading con- 
tribution to p for r — > is the lowest-energy part 1^^ . Its 
T-dependence cancels the T^ prefactor in Eo. 1(78(1 . and 



p (X JT — const. 



(83) 



Hence, the parameter p comes out as a (dimension- 
less) finite constant. With p in Eo. 1(66(1 being constant, 
the dynamical susceptibility l(65a(l immediately assumes 

with 

: UJ/UJO 

(84) 



a scaling form, x(q, u) = ^ [l + x"^ — py^ — iy\ 
length and energy scales ^ and ujq in x — , y 
given by 



Wo oc T/^ , C cx exp(0.63 J/T) 

The former has been given in Eo. 1(70(1 . the latter will be 
derived below. 

The structure factor defined in Eg. ((51(1 then takes the 
scaling form Eqs.(|3J), (@J), with the scaling functions 



ip{x) 



1 



1 + x^ 



, '^{x.y) 



l + x^ 



(1 + — py^)^ + y^ 



(85) 



The precise value of p is not known from the above 
calculation. The 'shape function' $(a;,y) is plotted in 
the inset of Fig. |51 for p = 1/4. It consists of a single 
peak at cij = for wave vectors in the relaxation regime 
X = <^ I , which splits into two peaks located sym- 
metrically to cj = for a; = fc^ > 5: in the regime of prop- 
agating modes. The boundary between these regimes is 
given by 9^$(x, y) |^_g = and reads x — ^(1 — 2p)/2p . 
The scaling function $ compares well to the numerical 
result shown in the main Figure|51. Apparently the choice 
p = l/A corresponding to a; = 1 is consistent with the nu- 
merics. For energies uj > too the damping (i.e., the width 
of the peaks for fc^ > 1) comes out larger in the numerical 
calculation, supposedly due to higher-order terms omit- 
ted in the ansatz ((66(1 for i?(w) in Ea. ((65all . The 'spin- 
wave' dispersion y{x) = uj{x)/ujQ of the peak-maximum 
is determined from dy^{x,y) — . At large wave vectors 
X it becomes linear, Q/ujo — ^^1 \fv- 
Finally the correlation length already stated in 1(84(1 has 
to be derived self-consistently from Eq. ( |65g| ) and (p!i^ . 
The integral in Eq.| |65g| is split into three parts, with 
the tanh(£/2T) approximated appropriately. 



H(0) = 



2^ (•QWo 

4t/ 



1 r 

piow{e) + j?f; depksie 
Phi{s) 



de 

2 Jt S 

H(0)^„ +H(0)T + n(0)j 
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With piou! from Ea. H74() inserted, the low-energy part 

becomes U{0)u;o = jijT^ !o '^VPiv) « i^o/T)'^/J « 
l/(^^ J) . Its contribution to the self-consistency condi- 
tion ((23 can be ignored, since 11(0) = 1/4 J + 0(C"^) . 
In n(0)T we can again use Ea. l(5n|) . i.e., phi{s) = , 
leading to n(0)T oc JT'^/u)'^ - T'^ / , which also can be 
ignored at small T <^ J . The dominant contribution to 
n(0) comes from high energies: Using p^i , Ea. H77|) . for 
LO > T results in 

1 r'^de 1 
2Tr J rp e uj^/e-e-ijJ 



5t^J J4T/3J {g'^-u^Y + u^ 

where g'^ = {4,Ld/iJf = ^'^\n{T /ljq) . The above ex- 
pression has to fulfill n(0)j = 11(0) = 1/4 J, giving an 
implicit equation for g^ and therefore ^ oc TjujQ. The 
exponential part of ^ can be extracted using T ^ 
and wc — > 00 in the integral boundaries. We find 
g^ = 1.07, and from Ea. (|70|l the correlation length is 

^(T) = yXexp(a J/T) , a = 37rgVl6 = 0.63 . 

VI. QUALITY OF THE CONSTRAINT 

The representation of spin operators in auxiliary par- 
ticles, Eq.lO, is valid only if the Fock space of canon- 
ical fermions is restricted to the subspace with Qi = 
Set fiafia = 1 on cach lattice site i . We have argued 
in Sects. ^ lllll IIVI that the local charge density Qi is 
a conserved quantity in the approximations considered 
in this paper, drQiir) = 0. In particular there is no 
spontaneous breaking of the corresponding gauge sym- 
metry in the 1/2-filled Heisenberg model. However, we 
have replaced the local constraint Qi = 1 by a global one, 
{Qi) = I 'i^ Qtot — J2i^ Qi — , since in that case we 
were able to use standard many-body techniques. 

The true ground state of the model Hamiltonian be- 
longs to the physical sector of the Hilbert space and 
thus features a total charge Qtot = with a homo- 
geneous density Qi = 1 . With just the global constraint 
Qtot = Nl in effect, this ground state might be (almost) 
degenerate with multiplets of unphysical states where 
Qtot is distributed inhomogeneously on the lattice. In 
that case the unphysical states mix spectral weight into 
correlation functions like the spin susceptibility. How- 
ever, one can hope that the separation in energy of these 
unphysical states from the ground state is finite, and at 
T ^ the correlation functions actually measure phys- 
ical processes. Since the Qi are conserved, the excited 
states generated in the correlation function (e.g., through 
a spin flip) remain in the physical sector. 

The gap of the lowest unphysical state to the ground 
state will depend on the physics of the latter. For the 
AF ordered state of the mean-field theory (see Sect. Illlj) 
chances are good: The Neel state \ AF) = \ AF; {Qi = 1}) 



is compared to some state \AF') = \AF; Qi = 0, Q/' = 
2, {Qj = Ij* 7^ ^)''}) I where a fermion has been trans- 
fcrcd from some site I to another site /' . Thus the spins 
on sites I and I' have disappeared, and the unphysical 
state \AF') is higher in energy by 2J, equivalent to 8 
lost J-bonds. The gap 2 J is just the charge-transfer gap 
liniT^o 2|m^| = 2 J in the spectrum Ea. (|17|l of the aux- 
iliary fermions. 

The above reasoning can easily be verified by calcula- 
tion, at least for the mean-field theory. The local charge 
(constraint) fiuctuations on an arbitrary site i are 

{{AQif) = {Q,Q^)-{Q,){Q,) 

= - ^ E ^ E {lu; + iv) G, (zc.)] 

involving a charge-response bubble built similar to the 
spin response Ea. l|19|! ^". With the fermion Green's func- 
tion in mean field, Ea. (|15|) . this becomes 

((AQ,)') =2/(m^)[l-/(m^)] (86) 

for i G A or B sublattice. Above the mean-field tran- 
sition temperature, T > Tn ~ J/2 , we have tua — 
and therefore {{AQi)^) — 1/2. That is, the constraint is 
significantly violated. However, in the AF ordered phase 
T < we have \mA\ > and charge fiuctuations are 
suppressed; at T <C TV Ea. l(5Bjl leads to activated be- 
havior, {{AQi)^) = 2exp(— 2J/r), revealing a gap 2 J 
between the Neel ground state \AF) and the lowest un- 
physical one \AF') . At T = the constraint is exactly 
fulfilled. 

For the minimal self-consistent approximation we do 
not have an expression for the free energy^^ and thus do 
not know how to calculate ((AQ^)^) consistently. An al- 
ternative test for the constraint is the local spin moment, 
which reads in the paramagnetic phase (Sf ) = , 

Since our system of interacting auxiliary fermions repre- 
sents a pure spin-1/2 model, we have ((Sj)^) — S{S+1) — 
3/4, and the local static structure factor should be 
^loc = 1/4 • Fluctuations of the fermion occupancy Qi 
will reduce from this exact value. At high temper- 
ature T ^ J the susceptibility Ea. H18|) is given by the 
simple fermion bubble Ea. H19(l . i.e., 

s?oc = ^Enr(*^) = lif^M) ihifl) = \ (87) 

iu 

where Eq.^ and (^/.t) = (4 /a) = (Q«)/2 = 1/2 
has been used. That is, constraint fiuctuations reduce 
the spin moment to 50% of its physical value. Note that 
this holds for any density of states of the fermions, in- 
cluding the case of free spins (J = 0). When temperature 
is lowered below ~ J in the interacting system, we ex- 
pect the spin moment to be restored as a consequence 
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FIG. 11: Influence of constraint fluctuations on the local 
spin moment S"^^ = {{Sf)^) , calculated numerically using 
Eg. lBH . Shown is 45'foc BiS function of the correlation length 
^(T) for the MSCA (crosses) and the fully self-consistent ap- 
proximation (diamonds). For ^ > 10 the MSCA reaches 85% 
of the exact value 4Sfoc = |'S'(5' -I- 1) = 1 . The dashed line 
indicates the high-T limit 1 /2 . 



of a strong reduction of ((AQi)^) in the short-range or- 
derd regions, ^f^^ has been calculated numerically using 
Ea. (|BT|l . the result for the MSCA is shown in Fig. ITTl as 
crosses. In fact the spin moment saturates to ~ 0.85 . 
The fully self-consistent approximation, which does not 
show short-range order, reaches ~ 0.68 (shown as dia- 
monds in Fig. Ill|l . 

For comparison we also discuss the spin moment in 
mean- field theory: At high T > T^r holds, and 

we have a reduced moment as above, ((S^)^) = 3| . At 
T <C Tjv , in the presence of long-range AF order, the 
spin moment is 



st\\ 
loc 



st± 
loc 



1 



q 



In the longitudinal direction just the condensate con- 
tributes, S'fgj,' = {S\)'^ , since fluctuations have vanish- 
ing weight, {S^SDconn ^ exp(— J/T) . The transversal 
contribution Sfl^ measures the spectral weight of spin 
waves and reads, from Ea. (|28|l . 



1 



loc 



1(^1)1(1 + 2e) 



1 



1 + 2(1-H2e)] 



with \{S\)\ = 1/2 + 0{e--'/'^) and the e introduced in 
Eq.(|231l . At T = the loop integral e = 0.197 is finite in 
2D . The spin moment is fully restored, the contribution 
of (transversal) magnons is even overestimated in spin- 
wave theory. 



VII. SUMMARY AND CONCLUSION 

This paper presented a study of critical fluctuations 
in the quantum Heisenberg antiferromagnet in two di- 
mensions. Our starting point is to re-write the spin-1/2 
operators in (unphysical) fermionic auxiliary particles, 
which enables the use of standard many-body techniques 
based on Feynman diagrams. Although the Heisenberg 
and the <- J- model have frequently been studied on mean- 
field level, the problem of strong long-range fluctuations 
in the two-dimensional plane has not been considered be- 
fore by this type of approach. 

The ground state at zero temperature T = shows 
long-range antiferromagnetic order and can be treated in 
mean-field theory and perturbative (in 1/5) corrections. 
The auxiliary fermions form a simple spin-density-wave 
state, their density of states consists of two delta peaks 
separated by a large gap 2 J . Results for magnetiza- 
tion, spin-wave velocity and -weight have been derived 
in Sect, mil . We found even quantitative agreement with 
conventional spin- wave theory. 

At finite temperature the perturbation theory in 1/5* 
breaks down, and the susceptibility has to be treated 
self-consistently in order to suppress the mean-field tran- 
sition temperature Tat = J/2 down to zero, as required 
by the theorem of Mcrmin and Wagner. In Sect. CSI 
we considered the simplest skeleton-diagram expansion 
in the renormalized fermion-propagator that provides 
this feature. Within a static (classical) approxima- 
tion (see Sect. IIV A|l the self-consistency problem can 
be solved analytically, leading to a correlation length 
^(T) cx exp(a J /T) and an exponentially small energy 
scale Wo esc , in agreement with the existing theo- 

retical and experimental work. The so-called dynamical 
scaling relationii is also fulfilled. However, the dynami- 
cal exponent z = 2 is incorrect, and the spin-excitation 
spectrum lacks a spin-wave like propagating mode. The 
analytical results have been confirmed by numerical cal- 
culation in Sect. lIVBl . 

We argued that the lack of 'spin waves' is due to a 
large density of states p{uj ~ 0) ~ 1/J of the auxiliary 
particles at the Fermi energy w = , which causes an ex- 
treme overdamping. It turns out that this problem can 
be cured by reducing the amount of self-consistency: In 
the skeleton-diagram expansion both the susceptibility 
and the Green's function of fermions are subject to their 
respective self-consistency loops, while being coupled to 
each other through the fermion self-energy. In the MSCA 
presented in Sect.0, which is not based on skeleton di- 
agrams, the additional loop for the fermion propagator 
is avoided. The results obtained from this approxima- 
tion are very satisfactory. In particular we find f (T) cx 
^ exp(0.359 J/T) and the energy scale wq 1/C j i-^-, the 
dynamical exponent is z 1 . The spin-excitation spec- 
trum now features dynamical scaling as well as a prop- 
agating mode for fc^ > 1 . k denotes the offset from the 
Neel wave vector. The MSCA equations have also been 
solved by analytical calculation in Sect. IV CI . confirming 
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in particular the energy scale luq = •\/3/27rT/^ . In the 
MSCA the spectral function p{uj) of auxiliary fermions 
features a strong suppression (pseudo gap) at the Fermi 
energy. The spectral weight at a; = is exponentially 
small, p(0) oc l/J^, and at low energies p{uj) shows scal- 
ing, P{UJ) = jfj:p{uj/uJo) . 

Our MSCA is technically and physically related to 
earlier work on low dimensional superconductors and 
the 2D Hubbard model, where long-range fluctuations 
of the superconducting or magnetic order-parameter am- 
plitude cause a pseudo gap in the single-particle density 
of states. References have been given in Sect. lV Al . Sup- 
posedly there is no ^-functional available^ which could 
have served as a systematic basis for the approximation 
scheme. However, these schemes follow quite naturally 
from a sequence of equations of motion in combination 
with a cumulant decomposition of higher-order vertex 
functions. This type of method has been used in the 
above-mentioned studies on superconductors and is ap- 
plied to the present problem in App. 

During all calculations performed in this paper the 
auxiliary-particle constraint Qi — fjafia — 1 has 
been treated in a mean-field fashion, {Qi) = 1 . Al- 
though it turns out that the fermions remain strictly 
local at 1/2- filling (no intersite hopping), on-site fluctu- 
ations ((AQi)^) into unphysical sectors of Hilbert space 
are not a priori suppressed. We did not attempt to im- 
prove on the constraint, e.g., by including a fluctuating 
Lagrange multiplier Ai(T) in the action Eq.® via the 

term Jq dr J2i ~ Qi (''')] • However, as has been 
explored in Sect. IVII . the actual amount of constraint 
fluctuations at low temperature depends on the physics 
of the ground state. Measuring the local spin moment 
((Si)^) in the MSCA revealed reasonably small fluctua- 
tions {{AQiY) at T ^ J . This is due to the magnetic 
order appearing on short length scales. At T = , with 
long-range order present, we even flnd ((AQ^)^) = in 
mean-field theory. 

A major question remaining open is whether the the- 
ory violates spin-rotation symmetry. In Sect. IIV CI it 
has been shown that a straight-forward application of 
the conserving-approximation scheme in two dimensions 
leads to a spurious phase transition at some Tjv ^ J , 
because the susceptibility is modified by additional ver- 
tex corrections that are not part of the self-consistency 
loop. However, since approximations constructed from 
this scheme do respect conservation laws only in the sense 
of a sufficient condition, it is worth testing our MSCA for 
spin-rotation symmetry 'by hand', e.g., through confirm- 
ing the Goldstone peak at q = (tt, tt) and w = in the 
long-range ordered state at T = . This requires a self- 
consistent calculation at (S) ^ , which we postpone to 
future work. 

We think that the auxiliary-fermion approach proved 
quite useful for dealing with strong critical fluctuations 
in the two-dimensional quantum model. Using standard 
many-particle techniques a suitable self-consistent ap- 
proximation can be found and treated with moderate ef- 



fort, even by analytical means. The approach may serve 
as a starting point for studies of strong fluctuations in 
the J-model at small doping levels. This would require 
to add bosons representing the doped holes (standard 
'slave-boson' formulation). Also the connection of the 
present study to the more general problem of pseudo-gap 
formation from long-range fluctuations might be helpful 
for the development of theories for cuprate superconduc- 
tors. Other possible applications include spin systems 
like the 2D Heisenberg model with frustration (e.g., the 
Ji- J2-model) or different lattice geometry. 
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APPENDIX A: SELF-CONSISTENCY 
EQUATIONS 



1. Equations for Numerical Solution 

For a solution of the Eas. (|35|l by numerical iteration, 
retarded propagators are introduced by analytic contin- 
uation to the real axis. 



A{iv) ^ A{uj + iO+) = A{uj) = A'{uj) + iA"{uj) 



where A = x^^tD^Yi . In order to avoid the numerical 
treatment of Bose functions g(w) we introduce the struc- 
ture factors 



5(q,c.) = [l+g(a;)]x"(q,^) 
S\uj) = [l + <?(u;)]H"(u;) 
U{u:) = [l+g{u:)]D"{u:) 
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As will be shown further below, the Eas. H35|l can now be 
written as 



U"{u; 



U{uj 



|i + j(q)nHP 



= -Pi ds 



n"{s) 



- / dep+{e)p (e-uj) 

/oo 
de [U{e)p'{e + Lu)+ 

+Ui-e)p+{e + u;)] 



TT 



de- 



7-4,7 |i + en(tj) 



(Ala) 
(Alb) 
(Ale) 

(Aid) 

(Ale) 

(Alf) 

(Alg) 
(Alh) 



with the density of states for the nearest-neighbor inter- 
action J(q) in two dimensions, 



. (27r)2 



^(e - 2J[cos{qx) + cos{qy)]) 



The temperature enters through the Fermi function /(w) 
in p^ , p^ . The output from the numerical iteration of 
these equations is the dynamical spin-structure factor 
5(q, w), Ea. l)Ala|l . and the spectrum of the auxiliary 
fermions p(a;) = [p{uj)^ + p(w)^] , Ea. ljAlel) . Further- 
more, the static structure factor (equal-time correlation 
function) is obtained using 



1 



div S'(q, iv) 



(All) 



In order to save computation time, the wave-vector 
integral in Ea. ljAlh|) is done analytically: J\f is approxi- 
mated by its value at the band edges, A/'(e) ~ ^^9(4J — 
|£|) . That is, the log-singularity at e = is ignored, 
which is anyway suppressed by the factor in U . The 
integral is now easily done, with the result (not writing 
the w-arguments) 



U 



K 



50 



7r4jn- 



8 J 



■In 



(M 



(n 



ir\2 



[M-f + {Wf 



2{IV 



tl\2 



n 



K 



(Alj) 



[arctan(M+/n") - arctan(M-/n")]/n" 

n'±4j|nP 



At 1^1 —> 0, where 11" (w) cx a; ^ , we use K{ijj) 
[l/M-{uj) - l/M+{uj)] . 



The derivation of Eas. ljAl|) goes as follows. Using the 
spectral representation of the fermion propagator G{iuj) 
in Matsubara space. 



G{iuj) 



, lUJ — £ 

the fermion bubble 135al) becomes 



At the real axis, w — > a; + iO+ its imaginary part is 
multiplied by [l-|-(7(w)] to form the 'bare' structure factor, 

poo 

= "2 y ^dep(e)p(e-c.)[/(£)-/(£-a>)][l+5(c^)] 

Using [/(e) - /(e - c.)] [1 + g{Lo)] - - [1 - /(e)]/(e - u:) , 
Ea. (|Ald|l immediately follows. Ea. (|Alb(l is just a result 
of the definition of S^{uj) and time-reversal symmetry, 
x"(-q,-c^) = -X"(q,^), n"(-c^) = -IV" {uj) . Simi- 
larly the imaginary part of the self energy l)35d|l reads 



~E"(c.) 



deD"{e)p{e + u)[g{e) + f{e + uj)] 



With .9 + / - .g(l - /) + /(I + 5) and D"{e)g{e) = 
i:'"(-e)[l + .g(-e)], which follows from D"{^e) = 
-L'"(e), Ea. (|Alf|l is obtained. The real parts H' and 
E' are just Kramers-Kronig transformations. 



2. Equations for Numerical Solution: Minimal 
Approximation (MSCA) 

The set of equations representing the minimal self- 
consistent approximation is identical to Eas. (|Al|l . 
except that the convolutions in and E" are missing, 
i.e., Eas. (|Ald|l and (|Alf|) are to be replaced by 



= -[U{u) + U{-u)] 



(A2) 



APPENDIX B: LOCAL SPIN MOMENT 

The spin moment ((S^)^) = 3(5fS'f) = iSf*^ is mea- 
sured through the local static structure factor 



TT Al ^ 

q 



J-00 ^ J~4j |l + en(c^)|2 
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using Eg. ljAlall and the density of states for 2D, intro- 
duced in App. I A II . The latter can be written 



27r2j 



K{Tn)Q{AJ-\e\) , m=l-{e/Uf 



with the complete elliptic integral K{m) = dt[{l — 
t-^){l — mt^)]^^^^ . In contrast to the calculation of U{lj) 
in Add. I A ll above, it is not suitable to approximate A/'(e) 
by a constant, and Sf^^ is computed numerically from 
Ea.dBlll . 



APPENDIX C: DERIVATION OF EQS.lHsl) FOR 
THE MSCA 

The equations Ht)5|l for the analytical treatment of 
the minimal approximation in Sect. IVT^ is derived from 
Eas. H60|l as follows: For wave vectors fc<Cl,k = q— Q, 
we have J(q) ~ —2dJ + Jk^ , and the susceptibility 
can be written 



J 



/ 1 



-2d] + 



1 



VJn(o) ' ' \jn{uj) jn(o). 

Introducing the correlation length ^ as in Ea. (|39|l . rep- 
resenting the static and n(ct;) containing the dynamical 
behavior of x i 

/ 1 



- 2d 



V Jn(o) 

it follows Ea. (|65a|) . with 



n(Lj) = n(tj) - n(o) 



Rio;) 



1 



1 



J Vn(cj) n(o) 



= 2df 



2djn(cj) 

1 + 2dm{ijj) 



Here 11(0) — l/2dJ has been inserted, omitting terms of 
0(^~2) . For dispersing peaks (i.e., propagating modes) 
to be able to appear in Imx for ^fc > 1 we may assume 
that \R{uj)\ - ^2;,2 therefore |n(cj)| - fc^ < 1 . In 
d = 2 this leads to Ea. (|65b|l . 

The fermion propagator Ht)Ue|l reads 



G{iuj) 



its spectral function p{e) has been defined in Ea. (|6& 
Now the irreducible bubble 161) all turns into 



(CI) 



1 

^ J — oo 



tanh(£/2r) 



At iv — !■ cij + iO+ , the imaginary part becomes Imn(a;) = 
Imn(a;) — ^ tanh([x)/2T)/9(w) , which forms the imagi- 
nary part of Ea. (|65c|) . The expression for 11(0) given in 
Eq. ( |65g| ) follows immediately with v — ,\i the particle- 
hole symmetry p{e) — p{—£) is exploited. The real part 
of Eg. j^ comes from Ren(w) = Ren(w) - n(0) , 



Ren(w) 



de tanh(£/2r)p(£) 



1 



1 



using again p(e) = p{—s) . 

The renormalized interaction (i.e., the local suscepti- 
bility) is also written in spectral representation, 



1 r°" 

D{iv) = - / de 

J-oo 



lmD{£) 



£ ~ IV 



(C2) 



Inserting this into the self energy Ij60d|) and performing 
the Matsubara sum results in Ea. (|65c|l . 

Finally, D as defined in Ea. H60c|l can be written 



1 



Here only the contribution from wave vectors q close to 
the Neel vector Q = (tt, tt) is reproduced accurately. It 
is then consistent to let J(q)^ ~ ■^(Q)^ = {2dJY , which 
leads to Eq.JesJ if d = 2 . 



APPENDIX D: DERIVATION OF THE MSCA 
FROM THE EQUATIONS OF MOTION 

We follow the line of Ref. and consider the general 
Hamiltonian H = Hq + Hi for interacting fermions, 

i/0=^e.ctc. , Hi = l V:':cl4'CsCr (Dl) 



r collects all single-particle quantum numbers; for the 
Heisenberg model Eq.(^ it is r = (i, a) , and the interac- 
tion matrix element reads TA^,'^ = J^^ icr"'"cr^'''^ . The 
commutator of a fermion operator with the Hamiltonian 
is given by 

[Cr , Hq] = Sr Cr , [c^ , = ^ K^^ra 4^ '^''2 ^ri 

'•l,I'2,''2 

(D2) 

The propagator for n fermions is defined as 

G„(12 . . . n; l'2' . . . n') = (-l)"(r,CiC2 . . . cj^, . . . cj,) 

with the short-hand notation ci = (ri) , cj, = c^., (t{) 
etc. The equations of motion for the 1- and 2-particle 
Green's function are given by (see, e.g., Ref. Flh. 

(a,, +e,JG(l;l') ^ {Tr[cr„Hi]{n)cl,)-S{l,l') 

(9,, +e,JG2(12;l'2') = -(T, [c,, , i/i](ri) cs^cj,) 
-<5(1,1')G(2;2') + <5(1,2')G(2;1') 

with (5(1, 1') = 5{ti — T[)6r-i^^r[ ■ Inserting the commuta- 
tor (|D2|I and applying J dl G(3; 1) to both sides of the 
resulting equations, we obtain after re-labeling 1 — » 1 , 
3-.1, 

G(1;1') = G(1;1')+ (D3) 
di'dld2'd2 G(l; l') V{1'1; 2'2) G2(l-2; 1'2'^) 
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Here l± = (n ,Ti±0+),Jdl = An , F(l'l; 2'2) = 

V^^'j,^M(t{ — ti)5{t2 — T2)6{ti — T2) , and G denotes the 
bare {V = 0) Green's function, which fulfiUs 

(a,, +e,JG(l;l') = -<5(l,l') 

Similarly, 

G2(12; 1'2') = G(l; 1') G(2; 2') - G(l; 2') G(2; 1')+ (D4) 
+ / drdId2'd2G(l;r)l^(l'l;2'2)G3(l_22;l'2'2^) 



The Equations (|D3p and (|D4|I are exact. Self- 
consistent approximations for the 1-particle Green's func- 
tion G result from a decomposition of the 3-particle func- 
tion G3 in the r.h.s. of Ea. (|D4|) into lower order propaga- 
tors. For that we note the generating functional for the 
n-particle Green's function^", 



T>[c, c] e 



A = j dl Ci{dr,)ci + J^dTH{T)+ J dl (ciT^i +r7ici) 



G„(l...n; l'...n') = 



1 



Z 6rii5rji' . . . Sr]nSr]n' 



+ 



• < 



2' 2 2' 2 
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< • < 
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FIG. 12: Diagrammatic representation of Ea. l)D7^ . which has 
been iterated. All diagrams corresponding to vertex correc- 
tions have been omitted, see text. Thin and thick continuous 
lines denote Green's functions G and G , respectively. The 
dashed line is —V, fermion loops contain a — sign. The ex- 
change diagrams, obtained via 1' <-> 2' , are not shown. 



The corresponding connected Green's function is defined 



G^(l...n;l'...n') 



!;2ji 



HZ) 



SriiSrji' . ..Sr]„Sr]n' 



(D5) 



i)=0 



Applying Ea. llD5l) for 71 = 1,2 leads to the familiar results 
G^(l;l') = G(l;l') and 

G^(12;1'2') = G2(12;1'2')- 

-{G(1;1')G(2;2')-G(1;2')G(2;1')} 
whereas the expression for n ~ 3 reads 
G^(123; 1'2'3') = G3(123; 1'2'3')~ (D6) 

- {G(l; 1') G2(23; 2'3') -G(l; 2') G2(23; l'3') + 

^^^^^^^^^"v^^^^^^^^"^ ^^^^^^^^^~*v ^^^^^^^ 

(rpa) (rpacx) 

-I- G(l; 3') G2(23; 1'2') + all other contractions} + 
+ 2{G(1;1')G(2;2')G(3;3')- 

- G(l; 1') G(2; 3') G(3; 2') -I- aU other contractions} 



The l.h.s. of Ea. (|D6(l represents effective 3-particle inter- 
actions and is ignored^^*^ . From the r.h.s. we keep only 
G3 and the terms (rpa) and (rpaex), and Ea. (|D4p turns 
into a Bethe-Salpeter equation, 

G2(12; 1'2') = G(l; l') G(2; 2') - G(l; 2') G(2; l')+ 

di'dId2'd2G(l;r) V"(ri;2'2) x (D7) 

X {G(i_; 1') G2(22; 2'2;) - G(i_; 2') G2(22; 1'2;)} 



This equation can now be iterated, leading to an expan- 
sion of G2 in powers of V through the series of particle- 
hole bubble diagrams shown in Fig. 1121 . 

Note that the bubbles contain one bare G (thin line) 
and one full Green's function G (full line), which is a 
direct consequence of the combination G G appearing in 
the first line of Eq. (|D4|I and accordingly in HD7|) . The it- 
eration of Eq. ljD7|l additionally produces diagrams which 
either contribute to vertex corrections to the simple bub- 
bles, or form self-energy insertions of Hartree type. The 
former have been omitted in Fig. 1121 . The latter contain 
loops of the form / dl'dl V{2'2; I'l) g(i; l'^) cx crTr[cr] = 
, g — G,G , which are zero for the Heisenberg model 
without magnetic order. For the same reason the other 
terms (except (rpa) and (rpaex)) in Ea. ljD6|l are not con- 
sidered. 

Finally, the G2 shown in Fig. E| is inserted in 
Ea. (|D3l) . omitting again all vertex-correction diagrams 
that emerge from the contraction of G2 in (|D3|I . Eq. HD3() 

now takes the form of a Dyson's equation, 

G(1;1') = G(1;1')+ (D8) 
di'dlG(l;r)S(r;I)G(i;l') 



where the self-energy S is identified as the series of dia- 
grams obtained from 



S(l';l)= / dl'dl y(l'l; I'l) G(l;l') 



(D9a) 
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with the renormalized interaction Eas. ljD8|l and (jD9p represent exactly the minimal self- 

consistent approximation given in Fig. [7|. 

y(l'l;2'2) = F(l'l;2'2)+ (D9b) 
+ / drdId2'd2y(l'l;ri)G(I;2')G(2;r)t/(2'2;2'2) 
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